Motivation and Implementation of Euler's Method 

1. What Would One Use Numerical Methods For? 

The graphical methods described in the previous session give one a 
quick feel for how the solutions to a differential equation behave; they can 
also be very accurate at predicting long-term behaviour, e.g. in the presence 
of funnels. However, when the 'medium range' solutions must be known 
accurately, and the equation cannot be solved exactly, nimierical methods 
are usually the best option. 

Also, even when an equation can be solved with an exact formula, it still 
might not be straightforward to compute values of a solution. For example, 
the equation y' = y with initial condition i/(0) =1 can be solved exactly: 
y{x) = e^. The niunber e is the value y(l). But how do you find out that 
in fact e = 2.718282828459045.... ? Here, too, you would use some kind 
of numerical methods. For ODE's the simplest numerical method is called 
Euler's method. 

2. Review: The Tangent Line Approximation 

Consider a function y{x), and a point {a,y{a)) on its graph. Call Ta{x) 
the function describing the tangent line to the graph at this point. The slope 
of the tangent line is y' (a) and it satisfies the equation 

Tfl(fl + h) = y{a) +y'{a) ■ h 

Suppose we are given a differential equation y' = f{x,y), with initial 
condition y(a) = c. One simple way to approximate y near a is to use the 
tangent line at a: 

y{a + h) ^y{a)+y\a)-h 

For example, consider the differential equation y' — y, with initial con- 
dition y(0) = 1. The tangent line approximation gives: 

y(l)«y(0)+y'(0)xl = l + l=2 

Of course, we know y{x) = e^, so y(l) = e 2.718. So, our estimate 
of 2 is rather crude. Suppose instead we were trying to approximate i/e 
(= 1.649...). We get 



^/e = y(l/2) ^ y(0) + y'(0) x 0.5 = 1 + 0.5 = 1.5 
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a better estimate. Graphically, we can see the tangent line separates 
from the curve. The bigger the step h in the x direction the farther 
away the tangent line is from the curve. 

The point is that the tangent line approximation works best for small 
h. The idea behind Euler's method is to use a sequence of successive 
tangent line approximations, each of them with a fairly 'small' h. 



3. Euler's method 

We want to estimate the solution (integral curve) to 
y' = f{x,y) passing through (xo/i/o)- It is shown as a curve in the 
picture. We first choose a step size, denoted h. Starting at (xo,i/o) 
we approximate the integral curve over the interval [xq, xq + h] by 
the tangent line, which has slope f{xo,yo)- (This is the slope of the 
integral curve, since y' = f{x,y).) This takes us as far as the point 
{x\,y\), which is calculated by the equations (see the picture) 

xi = xo + h, yi=yo + hf{xo,yo) ■ 

Now we are at (xi, 1/1). We repeat the process, using as the new approx- 
imation to the integral curve the line segment having slope f{xx,yi). This 
takes us as far as the next point {x2,y2), where 

X2 = xi+h, y2=yi + hf{xi,yi) . 

We continue in the same way. The general formulas telling us how to get 
from the (n — l)-st point to the n-th point are 

Xn = Xn-l+h, yn=yn-l+hf{Xn-l,yn-l) ■ (1) 

In this way, we get an approximation to the integral curve consisting of 
line segments joining the points (xo,i/o), (xi,i/i), ... as shown in the figure 
below. 
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We will call the line segments "Euler struts", and their imion the Euler 
polygon. It is an approximation to the integral curve y = y{x). 

In doing a few steps of Euler's method by hand, as you are asked to do 
in some of the exercises to get a feel for the method, it's best to arrange the 
work systematically in a table. 

Example. For the IVP: y' = — y^, y{l) = 0, use Euler's method with 
step size .1 to find y(1.2). 

Solution. We use f{x,y) = x^ — y'^, h = .1, and (1) above to find Xn and 

yn- 

n Xn y„ f{Xn,yn) hf{Xn,yn) 

10 1 .1 

1 1.1 .1 1.20 .12 

2 1.2 .22 
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